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We calculate the Aslamazov-Larkin term of the conductivity in the presence of a magnetic field 
in ■ applied along the c axis from the time-dependent Ginzburg-Landau equation perturbatively using 

two approaches. In the first a uniform electric field is explicitly applied; in the second the Kubo 
formula is used to extract the linear response. The former yields a version of the flux-flow formula 
for the uniform afo-plane conductivity, a xx (\i = 0), that holds to all orders of perturbation theory. 
Obtaining the same result from the Kubo formula requires considerable cancellation of terms. We 
O ' also use the Kubo calculation to examine the nonlocal afe-plane conductivity, a xx (k ^ 0) (where the 

»' | , cancellations no longer occur), as well as the nonlocal c-axis conductivity cr zz (k 7^ 0), and look for 

i 1 the perturbative precursors of the growing viscous length scales. In addition, we consider the effects 

of weak disorder — both uncorrelated (point defects) and correlated (columnar and planar defects) . 
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I. INTRODUCTION 



^3 

O, 

O , Recent experiments on YBCO have exploited nonuniform current distributions to reveal the nonlocal nature of 
the conductivity in the vortex-liquid regime. The conductivity is nonlocal if the current at a site r is determined not 
only by E(r) the electric field at r but also by the fields at other sites E(r'). The conductivity is always nonlocal on 
some microscopic scale, but the term "nonlocal" is generally reserved for cases in which the length scale involved is 
much longer than some underlying microscopic one. A nonlocal conductivity implies a nonlocal resistivity, though 
the length scales on which each varies may differ. Nonlocality may result if the "integrity" of the vortex along the c 
axis is maintained (i.e. not much cutting and reconnecting) as then a current applied at the top of a sample will yield 
a response not only at the top but also at the bottom. The interactions between vortices may also lead to nonlocal 



■ effects within the ah plane. p, 
Q\ \ The nonlocal nature of the conductivity has been studied by Safar et ala They injected small currents into the top 
of a sample with a magnetic field applied along the z axis. They extracted the current in one case from a second 
contact along the top and in another case from the bottom. In each configuration a series of voltage differences were 
taken. The results were then analyzed assuming a local though anisotropic conductivity (o~ xx 7^ o~ zz ) — the so-called 
modified Montgomery analysis. Safar et alH found an effective value for the ratio o~ xx /a zz for each configuration and 
discovered a discrepancy of order ~ 10 5 between them. The inadequacy of the local analys^ is taken as evidence for 
a nonlocal conductivity. Similar evidence for nonlocal effects in BSCCO has been reportedB. _ _ 
q \ It should be pointed out that the nonlocal interpretation has been questioned both on theoretical] and experimental!!! 
grounds. It should also be noted that the large effects seen by Safar et aL l i were in heavily twinned samples. Evidence 
for nonlocal effects in pure samples is somewhat weaker. Lopez et aln made a direct comparison of twinned and 
untwinned samples. In the first of the current configurations described above (current into and out of the top), they 
took voltage measurements across the top and bottom of their samples as a function of temperature. In twinned 
samples of YBCO, the resulting curves Vt op and V& t became almost indistinguishable at temperatures above the 
melting transition, indicative of a long c-axis length scale in the vortex liquid regime. In untwinned crystals, on 
the other hand, the curves only met at the melting transition. This result demonstrates the enhancement of the 
vortex-line integrity and hence nonlocality by the correlated defect. 1-. 

On the theoretical side nonlocal conductivity has been studied within the hydrodynamic frameworkB. We are 
assuming a linear but nonlocal relationship between current and electric field, i.e. an integral version of Ohm's law 

i^(r) = J a^(r,r')E v (r')dr', (1) 
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where o"^„(r, r') 7^ <7^S(r — r'). If the system is translationally invariant, Fourier transforming equation ([!]) results in 

j M (k) = «r^(k) £ w (k), (2) 

where nonlocality now implies cr^^k) ^ const. The emphasis is on the large-distance -behavior in real space and 
so presumably on the small-k behavior of the Fourier transform. Huse and Majumdaiij explored a hydrodynamic 
approach, in which the small-k expansion of er^^k) is truncated as follows: 

CT^k) = (7^(0) + Spapv k a kg. (3) 

In the example they work out in detail to explain features seen by Safar et ai, they include only one S: S xzzx . 
It models the "tilt viscosity" — which measures the influence of pancake vortices moving in one ab plane on those 
moving in another when they travel with different velocities.Q 

Stability requires that CT M ^(k) be a positive definite matrix. For the hydrodynamic form of a^{k), it implies, for 
example, that S xzzx and S xxxx mus t be positive. When these conditions are met, one finds that the resistivity decays 
with a length scale of order (y/S/a), which is in turn the length scale of the electric-field and current distributions .0 
(The hydrodynamic approach also predicts surface currents which supposedly should be spread out over some small 
length scale not accessible by that technique.) Recent calculationafm starting from the time-dependent Ginzburg- 
Landau (TDGL) equation find that many of the S"s do not have the requisite sign — at least not for clean samples at 
high temperatures. Those calculations were taken only to Gaussian order (lowest order in perturbation theory); the 
purpose of the present paper is to go beyond Gaussian order and also to explore the effect of various sorts of defects: 
point, columnar and planar. Recent simulations of the TDGL equations in two dimensions did find a region above 
the melting transition in which S xxxx > 0.EJ (The problem of calculating tip, voltage distribution whether or not the 
conductivities are of the "hydrodynamic" form will be presented elsewhere. Hj) 

The first order of business is to calculate the uniform conductivity a^(0). The resistivity of clean samples of 
YBCO in the vortex-liquid regime varies smoothly with temperature until it experiences a sudden drop generally 
associated with the first-order melting transition to a-pinned vortex solid. The temperature and field dependences 
in this smooth region agree with the flux-flow formulate which was originally conceived of in the vortex-solid phase 
where one pictures the translation of the entire unpinned vortex lattice. So how is it that the same physics applies 
in the vortex-liquid regime? The important feature of the flux-flow formula besides the absence of pinning by defects 
would seem to be its insensitivity to viscous effects. As it is the entire vortex system that moves, the amount of 
entanglement along the c-axis or the extent of crystalline order in the ab plane is irrelevant. It takes a nonuniform 
current distribution to probe such viscous effects. 

We will determine cr Mt ,(0) in two ways. First, we will consider the TDGL equation with an explicit uniform electric 
field applied, calculating the current and extracting cr^(0) directly. Second, we will calculate the conductivity from 
the Kubo formula, a version of the fluctuation-dissipation theorem which yields the linear response of a system not too 
far from equilibrium. The electric field does not explicitly appear in the latter. We will see that a delicate cancellation 
of terms arising in the Kubo calculation is necessary for the two approaches to agree, and we will demonstrate this 
cancellation to second order in perturbation theory for the conductivity of two-dimensional samples. We will then 
proceed to examine the nonlocal effects within the Kubo formalism, see how the k = cancellation breaks down at 
k ^ and show how viscous effects may enter. 

Next we consider the effect of disorder — first on the uniform conductivity. It would seem to have two competing 
effects. On one hand, defects pp. the vortices and thus might lower the resistivity. On the other hand, they disrupt 
the formation of a vortex latticelij which is more readily pinned than individual vortices, and thus they might raise the 
resistivity. In fact, both effects are seen in the experiments of Fendrich et alH3 in which point defects were induced by 
electron irradiation. Introducing the defects eliminated the sudden drop in resistivity at T c . At temperatures above 
T c , the irradiated sample had a lower resistivity; while at temperatures below T c , it had a higher resistivity. Another 
outcome of inducing the defects is that the temperature and field dependence is no longer of the flux-flow type. At 
low temperatures an activated form might be expectedEj, but such considerations are beyond the scope of the present 
work which employs standard perturbation theory. We will examine the effect of these defects in the weak disorder 
limit. 

Even the Jawest-order calculation of cr^^k, u>), the wave- vector and frequency-dependent conductivity, is 
cumbersome.B'EI To simplify the results we will concentrate on the leading behavior. In the presence of a mag- 
netic field, Landau levels provide a natural basis for calculating and describing various phenomena. Two energy 
scales naturally arise: the first is an, the energy of states in the n — or lowest Landau level (LLL); the second is 
the energy spacing between Landau levels given by Tiloq with loq the cyclotron frequency. We will gain an immense 
simplification by focusing on the regime in which an <C hu>o, the so-called LLL approximation. Arguments based 
on the renormalized values .of an and Tlloq suggest that such a separation of scales holds over a significant portion 
of the vortex-liquid regime.E^I Note that an <C Tiojq limit does not imply that only n = states are employed in a 
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calculation, although it usually does imply that the number of n ^ states is kept to a minimum. Considering only 
fluctuations in the LLL states serves not only to simplify the calculations but also to regularize them. 

Before presenting the calculations in detail let us highlight a couple key features. The Kubo formula for the 
conductivity involves the product of two current densities, and each current density involves a product of a \& and a 
where ^ is the superconducting order parameter. (See eqs. d34| ) and ( |35| ) below.) Thus, calculated in this way the 
conductivity is a "four-point" object, that is, it involves four iff fields. If these 'J's are expanded in the Landau-level 
basis, a crucial question arises: how many of these four 'f's can be in the LLL? Already we find a notable distinction 
between the c-axis conductivity, cr Z2 (0), for which the answer is all four and the a6-plane conductivity, <j xx (0), for 
which the answer is only two. As a result, in the LLL regime the characteristic time scale inherent in a zz is much 
larger than that in a xx . Another outcome is that in terms of LLL states (cr xx (0)) is a two-point quantity and as such 
is independent of viscous effects which are related to four-point quantities. I— ■ 

Understanding the effect of random pinning in the vortex-liquid regime is no trivial matter .113 When point disorder 
is added, translational invariance is destroyed. Moreover, while the division of cr xx (0) into n — and n > 1 parts 
persists, the n = component now becomes a four-point object sensitive to viscous effects. As a result the associated 
time and length scales may grow. The separation of n = and n > 1 states also plays a role in the distinction between 
c-axis length scales in materials with correlated and uncorrelated disorder. To each energy level there is associated a 
c-axis length scale, and the one linked to the n = level is much longer. It turns out that the c-axis length scale of 
a xx is controlled by the n = states in the presence of columnar and planar defects but by the n > 1 states in the 
presence of point defects. 

The rest of the paper is organized as follows. In the next section we lay the groundwork for the calculations that 
follow by constructing the Green's function and correlation function from the TDGL equation. In the section following 
that we calculate the uniform conductivity, cr xx (0), first by explicitly applying the electric field and then via the Kubo 
formula. We then move on to consider the nonlocal conductivity a xx ()s.). After that we examine the conductivity of 
films. We then consider the effect of weak disorder and finally summarize. 



II. TIME-DEPENDENT GINZBURG-LANDAU THEORY 



We begin with the usual Ginzburg-Landau free-energy functional J-[^f] for the superconducting order parameter 
^(r, t) in a magnetic field 



E 



(4) 



where A is the vector potential and Pj is the j component of the momentum operator given by 

i dr. 



(5) 



with e* = 2e. This combination arises for reasons of gauge invariance. We allow for two masses m x ,y — m a b and 
m z = m c to reflect to some extent the anisotropy of high T c materials; of course, the calculations can be extended to 
the Lawrence-Doniach model in which the layering is treated more explicitly. 

Next, we choose a simple relaxational dynamics for '5 given by the time-dependent Ginzburg-Landau equation 



1 fd ie*$(r,t) 
f\dt + n 



*(r,t) 



&F[<E] 



+ r)(r,t). 



(6) 



The field $(r, t) is related to the chemical potential which accompanies the time derivative in order to maintain gauge 
invariance; in what follows we will use the approximation that <&(r, t) is the scalar electric potential.li2rE3 The thermal 
fluctuations are represented by the noise ?7(r, t), which has zero average and S- function correlations 



( V *(r,t) V (r',t')) = ^8(r-r') S(t-t<). 



(7) 



The noise strength is chosen so that in the absence of a driving electric field, the distribution of ^'s, P^], evolves 
toward its equilibrium solution 



7>[*]cxexp j-^-L J d 3 r ^[*]| . 



(8) 
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We take T, the kinetic coefficient, to be real. Note that a complex kinetic-coefficient (F 1 — > F 1 + iA 1 in eq 



but T 1 — > T 1 in eq. (0)) is required to model the Hall conductivity er^.E] 
The TDGL equation is supplemented by 



V x V x A = /j J (n) + J (s) 
where J' n ' is the normal current given by 

j(n) = a (n) _ V(J) 
and j( s ) is the superconducting current given by 



dA 
~dt 



J\ s) ^t) 



2m, 



[Pij+Py] **(r 2 ,t)*(n,f) 



(9) 



(10) 



(11) 



with Pj is the momentum operator (|^) . Eq. (|^) is simply Maxwell's equation sans the Maxwell displacement current 
which is presumed negligibly small. Actually in this paper we concentrate solely on the TDGL equation, so that 
certain effects that enter with the Maxwell equations, such as backflowEjcil, will be missing. 

Let us consider an electric field applied along the x direction: $(r, t) = —E x and a magnetic field applied in the z 
direction: (A = Bxy). The TDGL equation then becomes 



Id 



i>(r,t) =77(r,t)-/3|V>(r,i)| 2 </>(r,t), 



where TL is 



H 



d 2 ( d ie*Bx 



h 2 8 2 
2m c dz 2 



2m a b \ dx 2 \dy fi 
For the calculations that follow we will need the Green's function which satisfies: 

G(r,t;r',t') = S(r-r') 8{t-t'). 



ie*Ex 

Th ■ 



(12) 



(13) 



(14) 



The Green's function serves as the inverse of the operator (r dt + Tl), allowing us to rewrite the TDGL equation as 
*(r, t) = J dr' J dt> G(r, f; r', t') [r?(r', f ) - /3|*(r', t')\H(r' , t 1 )] , (15) 

which we can write in a more symbolic form 

*i = Gi,arj2-j8Gi, 2 *2*a*2. (16) 
and solve — at least formally — by iteration 

*x = G li2 77 2 - pG lt2 G* 2}3 G 2A G2, 5 v* 3 V4V5 + 0(/3 2 ). (17) 

This expansion forms the basis for the standard perturbation theory in f3. 

We can construct G(r, t; r',t') from the eigenstates 4> n (k y , k z ;r) and eigenvalues E n (k y , k z ) of Ti., which are 



tfrnikyi k z ;T) — e 



e ik v y+ik z z 



n ( „ k y £ 



ivm ab t 
h 2 T 



and 



_ ,, . , h 2 k 2 fiu m ab v 2 ikyV 

E n(k v ,k z ) = ^—+a+ — + w ^+ n hu -—, 



(18) 



(19) 
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where 



e*B 



1/2 



LU 



e*B 
m ab 



and 



E 

V= B- 



(20) 



{I is the magnetic length and is associated with the distance between vortices; ujq is the cyclotron frequency; and v is 
roughly speaking the speed at which the flux lines move.) The functions u n (s) are given by 



u n (s) 



H n (s) exp{-s 2 /2} 



(£^2 n n\) 



1/2 



(21) 



where H n (s) are the Hermite polynomials. 

Note that in the absence of an electric field (v = 0) the energy does not depend on k y giving the characteristic 
large degeneracy of the Landau levels. Moreover, we then find in expression (|lj) the two energy scales mentioned in 
the introduction: 1) an = a + %uiq/2, the energy in the LLL at k z = and 2) hujo the energy spacing between levels. 
Since the temperature at which an — is where one expects the LLL modes to go critical within mean-field theory, 
it is taken to define the mean-field H C 2(T) line. 

Given the presence of the k z term in the energy, it is convenient to construct two c-axis length scales, one corre- 
sponding to each energy scale; they are 



2m c a H 



1/2 



and 



m c LO 



1/2 



(22) 



The former is the standard mean-field c-axis correlation length, which is temperature dependent; while the latter is 
a c-axis version of the magnetic length 



4 = 



m ab 



1/2 



(23) 



which is temperature independent. 

Prom the eigenstates and eigenvalues above we construct the following Green's function 

G(r,t;r',f)=r /^/ff / ^ exp {ik y (y - y>) + ik z (z - z') - iu>(t - t')} 



ca u 



kjy£ 'hl^fj^ 



i n ky£ zki>£ 



n=0 



TE n (k y ,k z ) - ilu 



(24) 



where k v — m a bv/h 2 T. One might notice that when v ^ the arguments of the u's are not complex conjugates with 
x <->■ x'; that is because the operator H is not Hermitian. The real issue is that eq. ( pT| ) is satisfied. 

Various representations of G(r, i; r', t') are useful depending on the calculation in question. Another useful form of 
G(r, t; r', t') can be derived by performing the integral over ui and using the identity 



n=0 



H n (x)H n (y) 
2 n n\ 



,l2xyt-(x 2 +y 2 )t 2 ]/(l-t 2 ) 
(l-i 2 )V2 ' 



(25) 



ft leads to 



G(r,t;r',i') 



4tt£ 2 



exp ■ 



dk z 
2n 

-coth 



Thuj (t - t') 



ii -l 



sinh 



Thuj (t - t') 



exp ^ -r 



2 1,2 



2m r 



2^2 



2m a b 



[{x-x'f + {y-y' + v{t-t')f 



U 2 



(*-*') 



exp ■ 



i(x + x' - 2ikJ 2 )(y - y' + v(t - t')) 
2P 



+ ik z {z- z') } G(t-t'). 



(26) 
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In this form there is no longer any summation and the spatial dependence is simply that of a Gaussian, but the price 
paid is in the hyperbolic time dependence. One can see in the combination (y — y' + v(t — t')) the tendency for the 
vortices to move in the negative y direction under the influence of the Lorentz force. 

After the Green's function, the next quantity we will need is the correlation function C(r, t;r',i') = 
(^(r, t) x S*(r',t')). There are two complications here: the operator H is neither Hermitian nor translationally in- 
variant. As a result some of the formulas relating C and G that we have become accustomed to are inappropriate; it 
is best to resort to the definitions. Substituting the expressions for ^(r, t) and ^*(r',t') given by eq. ( |l5| ) into the 
definition of C(r, t; r', t') and performing the noise average yields 

C(T,t;v',t') = / dt " ( dr" G(r,ty,t")G*(r',t>;T",t") + 0(P). (27) 



Many of the calculations that follow will use the Kubo formula in which there is no explicit electric field (v = 0). 
Then the expression for the correlator simplifies somewhat, becoming 

a=o(r,t;r',*') = 2fc B Tr J g J ^ / ^ e xp {ik y (y - y') + ik z (z - z') - iu(t - t')} 



^ r 2 £2(fc 2 ) + ^ 2 • (28) 



U 

X 

n=0 



Now let us move on to using these expressions to calculate the conductivity. 



III. THE FLUCTUATION CONDUCTIVITY 



Having the Green's function and correlation function, we are now ready to calculate the conductivity. We will focus 
on the conductivity due to fluctuations in the superconducting order parameter — the so-called Aslamazov-Larkin 
term. The normal contribution <j( n ' must be included separately, and we are neglecting other possible contributions 
such as the Maki-Thompson or density-of-states terms.E3 

Recall that conductivities take on a form 

e* 2 Af t 

a<x — , 29 

m* 

where e* and m* are the effective charge and mass respectively, r is a characteristic time scale, and Af is the carrier 
density. Since we are interested in the conductivity due to superconducting fluctuations (l^l 2 ) will serve as the carrier 
density. An examination of the TDGL equation reveals that the kinetic coefficient T has dimensions (energy xtime)^ 1 ; 
so that time = fx energy)^ 1 . The question then becomes: what is the appropriate energy scale? Thus far two 
candidates — au and Tiluq — have emerged. 



A. The direct approach 



Now let us derive the flux-flow form of the uniform conductivity from a calculation with an explicit electric field 
applied. This approach has several advantages over the Kubo calculation. It extracts the conductivity from the 
current, a two-point object, instead of the current-current correlation function, a four-point object. Related to this 
point is the fact that many of the diagrams occurring in the Kubo calculation make canceling contributions. In 
addition, in the direct approach the a/tuvo — > limit (the LLL approximation) requires only LLL states, which is 
not true for the Kubo case. It can also be extended to include the nonlinear effects. (The advantage of the Kubo 
formalism is that it can be used to calculate the nonlocal conductivities (^(k) and cr Z 2(k).) 

From equation ( pT[ ) we see that the average current is obtained from the correlation function as follows 

2im a b \ox ox' J 

We will drop the superscript (s) hereafter. At the lowest order in perturbation in theory, we can insert expression 
@ into @) to obtain 



(30) 



G 



^ = k S^{l-i)I dt "I dr " G ^"^ «v,*w) 



(31) 



(r,t)=(r',i') 



We have verified a posteriori that in the direct approach the a/htoo — > limit is equivalent to using only n = states 
from the start. Restricted to n — states, the Green's function becomes 

F f dk r dk 

G (r 7 t;r',t') = ^= I j ^ J e W {ik y (y - y') + ik z (z - z>) - TE (k y , k z )(t - f)} 



xexp < 



(? 



e(t-t r ). 



(32) 



Notice that the derivative with respect to x in eq. ( j3l"l) pulls down a factor of —{x/l 2 — k y i — ik v ), while that 
with respect to x' pulls down —(x'jl 2 — k y 2 + ik v ). The wave vectors k y i and fcj,2 will eventually prove to be equal 
(momentum conservation), so when r is set equal to r' we find 



(■4) 



he*k v 
m a b 



C (r,t;r,t) 



m ab ThuJo 



l*o| 



(33) 



where Co(r,t;r,t) — (l^ol 2 ) denotes the order-parameter fluctuations in the LLL. (Strictly speaking it is (l^ol 2 ) i n 
the presence of the electric field, which is where the nonlinear effects would enter). Note that the conductivity has 
the expected form with the characteristic time given by r = 1/ThujQ and that any temperature dependence enters 
only through the carrier density and not through the characteristic time. 

So far this result is for the lowest order in perturbation theory; however, all terms in the perturbative expansion 
for C(r, t; r', t') begin with some G(r, t; ri, ti) and end with some G*(r', t'\ r n , t n ). It follows that the effect of taking 
the derivatives in eq. ( |30| ) will always be the same as it was above, and consequently, the flux- flow result, eq. (|3"3|), 
holds to all orders of perturbation theory in the LLL approximation. This result is remarkable in the simplicity of 
the relation between cr xx (0), a dynamic quantity, and (l^ol 2 ); a static quantity. One always expects the conductivity 
to be proportional to carrier density, but in this case the proportionality constant r turns out to be rather trivial — 
having no temperature dependence. 



B. The Kubo formalism 



We now turn to the Kubo approach. For a system not too far from equilibrium the conductivity can be calculated 
from the Kubo formula 

^(k,w) = 2j£r/ rf ( r -0 / e ik -( r - r ')--(*-*') (jW(r,t) J^(r',t')) c , (34) 

where the suffix c denotes the "connected" piece and where this expression allows for a frequency dependence. The 
Kubo formula relates the conductivity, a dissipative quantity, to fluctuations in an associated quantity, here the 
superconducting current. As such, it is a version of a fluctuation-dissipation theorem. It calculates the linear response 
to an electric field without explicitly applying one; the fluctuations above are those in the absence of an electric field. 
Inserting the expression for the current ( |il| ) into the Kubo formula (Q) yields 

^(k,c) = — f- f d(r-r') f d(t-t')e ik -( r - r ')-M*-*') 

8k B -l vn^rriu J J 



(Pi» + P^)(Ps» + Pl) (*(r 1 ,t)**(r a ,t)*(r 8 ,t')**(r4,0>« 



(35) 



r 1= r 2 = 
r 3= r 4 = 



One can see that calculated in this way the conductivity is a four-point object as mentioned in the introduction. 
Henceforth, we will consider DC results (u — 0) only and will drop the frequency dependence from our expressions. 

The aim now is to calculate averages of the sort (^(rj, tx)^*(r2, £2)^0*3, *3) v I / *(r4, t^)) c . At the lowest order in 
perturbation theory (the Gaussian term); one simply applies Wick's theorem 
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(*(n, t)¥*(r a , t)*(r Bj 0**('4, 0)c = (*(ri, t)**(r 4 , 0) (*(r 3j 0**(r a , t)), 



(36) 



retaining only the connected piece. The Gaussian term is represented diagrammatically in Figure An arrow 
corresponds to a Green's function, and a circle corresponds to the noise average. The combination arrow-circle- 
arrow constitutes a correlation function. Hence the conductivity at Gaussian order requires two correlation functions. 
Not drawn but of crucial importance are the momentum operators acting at sites (r, t) and (r',f') that make it a 
current-current correlator instead of a density-density correlator. 




FIG. 1. Gaussian order diagram for the conductivity. An arrow corresponds to a Green's function, and a circle corresponds 
to the noise average. Not drawn but of crucial importance are the momentum operators acting at sites (r, t) and (r',t'). 

Let us focus our attention for now on the uniform dc conductivities a xx (0) and a zz (0). Calculated in the LLL 
approximation and to Gaussian order they are 

4 ?(0) . £W> ( M (3 r) 

m ab \Thio a J 

and 

tt .dM(') (38) 
m c \8la H / 

We see here the characteristic time of the a6-plane conductivity is r a \, = l/Thuja while that corresponding to the 
c-axis conductivity is r c = l/STajj- 

Consider now how the Kubo formula reproduced the flux- flow result at Gaussian order. First of all the momentum 
operators intrinsic to a xx (0) act like creation or annihilation operators raising or lowering the Landau level with the 
result that one of the \P's at (r,i) is in a higher level; the same thing happening at (r',t') as well. Consequently, in 
the LLL limit, we have one n = correlator and one n = 1 correlator. With this in mind, we now look at the time 
integral in the Kubo formula, which becomes 

d(t-t') exp{-r[^(fc, 1 )+^(fc^)]|t-t / |}, (39) 

where k z \ and k Z 2 are the z components of momentum running through the n — and n = 1 channels, respectively. 
(Actually k z i = k Z 2 by momentum conservation.) Because Eo(k z %) is much smaller than Ei(k Z 2) we can drop Eo(k z i) 
from the integrand above. This last step is equivalent to using a static (equilibrium) n = correlator and hence 
produces a relation between the dynamic conductivity and the static density. In fact, the integral yields simply a 
factor 1 /Thioo — the characteristic time in the flux- flow formula. The flux- flow formula obtained in the direct approach 
above suggests that the Gaussian density (|\I/o| 2 ) should be replaced by the fully renormalized density (l^ol 2 ) but that 
the characteristic time T a \, = 1/TThjq should remain as is. The way to produce this result within the Kubo formalism 
is to renormalizc the n = correlator while leaving the n — 1 correlator essentially untouched. If this is true then 
diagrams that dress up or otherwise disrupt the bare n = 1 correlator should have no overall effect. 



C. The Hartree-Fock Approximation 

Our first attack on going beyond Gaussian order, i.e. incorporating the nonlinear terms, will be the Hartree-Fock 
approximation. At Gaussian order the superconducting fluctuations are 

(l*o| 2 ) = , kBT , f2 , (40) 
47r a H £ c t z 
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which would seem to imply that the conductivity diverges at mean- field, Hy-f (an = 0). This seeming divergence is 
eliminated when one adopts the Hartree or Hartree-Fock approximation.EifEj 

From the perturbation theory for W ( p"7| ) we can obtain the series for (\&(ri, t)^*(r2, t^fts, t')^*(r4, t')) and in 
turn the series for the conductivity. The new diagrammatic feature is a "quartic" vertex at which four lines meet 
with two arrows pointing in and two pointing out. Two diagrams occurring at the first order in are shown in Figure 
||. Because of the structure of the conductivity especially the momentum operators acting at (r,i) and (r',t'), the 
diagram shown in Fig. |^(b) has zero contribution at k = for both a xx and a zz . 





(a) (b) 

FIG. 2. Diagrams with one quartic vertex. Diagram ^(a) is included in the Hartree-Fock resummation. Diagram ^(b) has 
zero contribution at k = and is "down" in the calculation of a xx (\i ^ 0). 

Diagrams like that in Fig. ^(a) can be resummed by replacing an with a, where a is defined self-consistently as 
a = an + 2/?(|^ | 2 ), where what we mean here by (l^ol 2 ) is (l^ol 2 ) ( e q- (rhO|) with an — > a; all of which leads to 



(3 k B T 

a = a H + - . 41) 

Replacing an with a constitutes a resummation, the Dyson equation for which is represented diagrammatically in 
Fig. |3| It corresponds to a renormalization of the Green's function. 




FIG. 3. Diagrammatic representation of the Hartree-Fock approximation. The double lines represent renormalized Green's 
functions; the single lines represent bare Green's functions. 

In analogy with many-body physics, this replacement is called the Hartree-Fock (HF) approximation as it takes 
into account both the direct and exchange terms. In addition to eliminating a large number of graphs from the 
perturbation theory; the HF resummation replaces an which can be zero or negative with a which cannot — curing 
the divergence problem mentioned above. The removal of the divergence is connected to the fact that a large magnetic 
field effectively reduces the dimension of the problem bjj— two.BZI Moreover, the HF approximation goes a long way 
toward achieving agreement with measured conductivities t3 This success suggests that the same philosophy (aj? — > a) 
should be adopted when considering where the LLL approximation is valid. In the HF approximation an <S while 
Hujq remains the same. So whereas it originally appeared that the LLL approximation was valid near the mean-field 
H C 2 (T) line where an is small, it would now appear to be valid where a is small. But as already noted a dojea not 
change sign, it grows small only as an grows large and negative, i.e. below the mean- field H C 2(T) line. IkedaO has 
investigated renormalization effects beyond HF, finding that they are merely refinements to the HF considerations, 
suggesting that the LLL approximation has a substantial region of validity within the vortex-liquid regime. 

If we replace an with a in the Gaussian conductivities and look in the limit of an *C we find that 

lim a x J0) = e *)^i . (42) 
«««o xxy ' 2/3hBT v ' 

Note that <r xx (0) at this level of approximation already shows many of the features seen in the experiments in the 
vortex-liquid regime: the resistivity is linear in B with a zero intercept and extrapolates to the normal-resistivity at 
the mean-field H C 2(T) line, as can be seen for instance in the Fendrich et al. data prior to irradiation.Ell 







D. Beyond the Hartree-Fock Approximation 



While the Hartree or HF approximation goes along way toward describing certain properties such as the specific 
heat, it is totally inadequate for examining other features, such as the extent of-erystalline ordering within the ab 
plane. A much more sophisticated approach, such as the Parquet resummationc3, is needed for that. We do not 
provide such a scheme here, instead we examine a few diagrams beyond the HF approximation. We will find a 
cancellation among most of these diagrams at k = suggesting that the Hartree approximation is already quite good 
for the uniform conductivity. The absence of this cancellation at k 7^ 0, on the other hand, suggests the need to go 
beyond Hartree theory when considering the nonlocal effects. 

As we are now considering a perturbation theory about the HF approximation, the series is no longer a power series 
in (3. Rather it is a power series in x where 



(3k B T 



(43) 



16tt£ 2 £ c S 2 ' 

the dimensionless parameter introduced by Ruggeri and ThoulessH The self-consistent equation above, eq. (f4l|), is 
now rather compactly written as 

a H = &(l-8x). (44) 

In the expansion of the uniform conductivity a xx (0) there are no contributions of 0(x). Recall that of the diagrams 
at order (3 (Fig. ^) there are those of the HF type which have been absorbed into a, and the others make no 
contribution at k = 0. Some diagrams contributing to the conductivity at 0(x 2 ) are shown in Figures |], ^ and |[ 




FIG. 4. Diagrams contributing to a xx at order x 2 . The bold lines indicate n > 1 Green's functions; the thin lines n = 
Green's functions. Diagrams (a) and (c) are Green's function renormalizing; while (b) and (d) are noise-renormalizing — the 
difference being in the placement of the circles (noise averages). Diagrams (a) and (b) contribute to the flux-flow formula, (c) 
is "down" and does not contribute to the LLL approximation, and (d) is canceled by another diagram at k = 0. 
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FIG. 5. Vertex renormalizing diagrams contributing to a xx at order x 2 . Diagrams (a) and (b) are "down;" diagram (c) is 
zero at k = 0; and (d) is canceled by the diagram in Fig. g(d). 





FIG. 6. More vertex renormalizing diagrams contributing to a xx at order x . These diagrams do not contribute at k = 0. 



There are a number of considerations in addition to the usual ones of the topology and degeneracy to keep in mind 
when enumerating and evaluating the diagrams that occur in the perturbation series for a xx . Because the conductivity 
is a dynamic quantity, there is the distinction between Green's functions and correlation functions to bear in mind. 
We make this differentiation by including a small circle representing the noise average in the middle of a correlation 
function. Diagrams may differ in the placement of circles; for example, compare Figs, ^(a) and (b) or Figs ||(b) and 
(c). In addition, because the order parameter is complex, the Green's function carries a direction, indicated by the 
arrow. Hence diagrams may be distinguished by the placement of arrows around the diagram; for instance, compare 
Figs. H(c) and (d). Also because the order parameter is expanded in Landau levels, the Landau level of each line is 
another consideration. We make this differentiation by putting the higher Landau-level Green's functions in bold. So 
diagrams may differ in Landau level structure; for example, compare Figs. |J(a) and (c) or Figs. ||(a) and (b). On the 
other hand, in the evaluation of <r zz (0) all of the Green's functions are in the LLL. 

Recall that at Gaussian order cr xx (0) has one correlation function and thus two Green's functions in the n = 1 
state. It turns out that at all orders in perturbation theory the diagrams contributing to the LLL approximation have 
a maximum of two higher-Landau-level Green's functions. So diagrams like that in Fig. ^(c) and those in Figs. ^ 
(a) and (b) are "down" by a factor of a/hcoo and are not included to this order of approximation. Of the remaining 
diagrams it has been argued that those which have a simple n — 1 correlator (e. g. Figs, ^(a) and (b)) contribute 
to the flux-flow result and that the others {e.g. Figs. [|(d), ||(c) and (d) and g(a) and (b)) must have a combined 
contribution of zero at k = 0. When we evaluated these diagrams, we indeed found that the diagram in Figs. |^(c), 
|(a) and |(b) gave no contribution at k = and that the contributions due to the diagrams in Figs. ||(d) and Figs. 
5(d) were equal and opposite. We will demonstrate this explicitly in the 2D case. 

The diagrams in Fig. ||(a) and (c) renormalize the Green's function, and those in Fig. ||(b) and (d) renormalize 
the noise. The diagrams in Figs. ^| and ^ renormalize the quartic vertex. Note that on the n = 1 side, the Green's 
function renormalizing diagram was down whereas the noise renormalizing diagram was not. A similar phenomenon 
occurs when we consider the effect of disorder. 



11 



IV. WAVE- VECTOR DEPENDENCE 



Now let us turn to the nonlocal, k-dependent, conductivity. In some instances we calculate the full k-dependcncc 
while in others we restrict our attention to the coefficients in a small-k expansion 



o- M „(k) = 07^(0) + Spapv k a kp + 0(k 4 ). 



(45) 



The coefficients S^ap^ and the associated length scales can be related to the vortex picture. Since the vortices move in 
the y direction when the electric field and current are in the x direction, S xxxx is associated with shearingl and S xyyx 
with compression. S xzzx is related to the "integrity" of vortices (Huse and Majumdan3 call it the "tilt viscosity"). 



A. The Gaussian results 



The small wave-vector expansion of the conductivity at Gaussian order (within the combined HF and LLL approx- 
imations) is 



<r^(k) = 4^(o) 



(G)t 



(G) 

where we have factored out <j xx (0). First note that S xxxx and S xzzx are negative, i.e. they have the sign that cannot 
be handled using the hydrodynamic approach. Furthermore, the length scales multiplying k 2 and k 2 are magnetic 
lengths, i.e. they have no temperature dependence. 

(c) 

It is easy to see why (calculationally) the c-axis length scale in the Gaussian calculation of a xx is t c . The z- 
dependence has no effect on the Landau level structure; therefore, it remains true that cr xx (k z ) is comprised of one 
n = and one n = 1 correlator. We can take advantage of the disparity in "masses" (a <C hujo) by insisting that the 
external momentum k z be sent through the n = 1 channel. (This choice does not affect the outcome, it just makes 
it more apparent.) The internal momentum integral is dominated by the smallest poles, which are those associated 
with a, and the terms involving k z and Huq are left essentially intact. The full k z dependence at Gaussian order is 
thus 



Pk 2 x m 2 io 2 fk 2 v 



64<S 2 



o(k 4 ) 



(46) 



'£>(*.) = *£>(0) [1+^/2] 



(47) 



where the structure is simply that of [Ei{k z )]~ 2 with a — ► 0. 

The Landau level structure of S xxxx is slightly different. This time one correlator is in the n = level while the 
other is in either the n = 1 or the n = 2 level. For an expansion in k^. each additional power of k 2 requires one higher 
Landau level. But the important point is that there is no contribution with two n = correlators. The k x dependence 
at Gaussian order is 



4?(M = ^(o) 

The combined k x and k z dependences are given by 



1 



-k 2 J 2 /2 



k 2 J 2 /2 



-kle 2 /2 

n=0 



(n + l)(k 2 J 2 /2) n 
n\ \{n + 1) + fc^2/2] 2 



(48) 



(49) 



As opposed to S xxxx and S xzzx , S xyyx is positive at Gaussian order. In fact if we had included the next term in 
the u/Tiusq expansion we would S66 tllctt S X yy X Cllcill ges sign as T is lowered.El As T decreases the associated length 
scale associated with S xyyx , £± cx Tiwo^/cif, increases. A T-dependent length scale here is somewhat surprising. In the 
vortex picture S xyyx appears to be related to compression, but one might expect that the compressibility to be pretty 
much the same for the liquid and the solid, and so relatively T- independent. The Gaussian calculation is in conflict 
with this expectation. If this length scale found in the transverse conductivity does indeed increase as T decreases, it 
will prove interesting to compare it to similar growing afo-planeJength scales, for example, the phase coherence length 
and the length osier which density-density fluctuations decayed, which recent Monte Carlo simulations suggest grow 
in the same wayEil. 
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In the scenario in which we apply current and wish to extract the voltage distribution, it is the characteristic 
lengths of the resistivity which are important. This is why the signs of the S's are so crucial. They determine the 
pole structure of the nonlocal resistivity p(k) and consequently the length scales of p(r). Consider, as an example, a 
c-axis conductivity of the form 



CTzz(k z ) 



= „(*») 



+ ai s J(0) + S zzzz k 2 z 



0(k A ), 



(50) 



where a local normal conductivity has been included. If S is positive, the hydrodynamic approach^] yields a length scale 
[S/ (a^ +cr^ l - ) )] 1 / 2 . If S is negative, a Pade approximant approacha produces a length scale [er^ 1 5* | / cr^ s ^ (crW +a^)] 1 / 2 
which is much smaller than the previous one because of the factor {a^ / a^) 1 / 2 . Returning to the case of S xyyx , 
if it were to change sign from positive to negative as more diagrams are included, the associated length scale would 
then be short which may be consistent with the expectations of the vortex picture. With such dramatic consequences 
regarding the nonlocal behavior of the resistivity, one might expect that changes in the signs of the S"s would have 
experimental signatures. 

The difference between S xyyx and the others is that all four ^'s in the calculation of S xyyx can be in the LLL. 
Applying the definition of S xyyx , namely 



S x 



ld 2 a xx (k) 



:,,,J ~ 2 dkl 



k=0 



to the Kubo formula (eq. (pq)) and restricting to LLL states yields 



Jxyyx 



16k B Tm 2 



J d(r - r') J d(t - (y - y') 2 [p lx + P 2 * x ) [p 3x + P* ix ) 



(* (n, t)*5(ra, t)*„(r 3 , t')%(r 4 , f)>« 



(51) 



(52) 



This expression involves only LLL states is hence accessible to simulations and other methods which use only LLL 
states.H 

Another quantity that involves only LLL states is the c-axis conductivity, a zz (k), at Gaussian order it is 



(G), 



rt k 2 L i 2 /2 



l-(l + fc^ 2 /4)-V 2 
kit 2 /* 



(53) 



where k\ — k 2 + k 2 . Since the calculation involves only n = states, the c-axis length scale is necessarily £ c . If this 
result is expanded one finds that none of S's has the sign required by the hydrodynamic approach. 



B. Beyond Gaussian 

As already mentioned the HF approximation is inadequate for determining the extent of crystal ordering within 
the ab plane and as such is also inadequate for examining shear effects. To see if S xxxx has a long length scale or 
changes sign at low T, we must look beyond the HF approximation. The same is true of S xzzx . 

We argued that the uniform conductivity, a xx (h — 0), could be obtained from the subset of diagrams that renor- 
malize the n — correlator and leave the n — 1 untouched beyond the HF resummation. In the LLL limit the fully 
renormalized n — correlator C n= o(r, t; r', t') has exactly the same dependence on x, x' , y and y' as its "bare" version; 
thus, the fc^-dependence of this particular resummation is exactly the same as that in the Gaussian calculation (eq. 
(fl8|)). Furthermore, the external momentum k z can still be sent through the n = 1 channel, yielding again the de- 
pendence seen in the Gaussian calculation (eq. (|47|)). Resumming this subset of diagrams only alters the magnitude 
&xx(0); the k x and k z dependence remain the same as they were in the Gaussian calculation, i.e. 

ag F \k x> k z )=a xx (0) 

where 



-kit: 



72 E; 

n=0 



(n + l)(fc^/2)" 
\[{n+l) + mi2f 



(54) 
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= ,55, 
m ab l hu) 

and where 

So for this particular subset of diagrams, we only have to do the perturbative expansion for the static quantity (|^o| 2 )- 
Consequently, if S xxxx and S xzzx are going to change sign at low T and give long length scales, then other diagrams, 
such as the vertex-renormalizing diagrams (see Fig. [5]), must be important. 

Recall that in the a <C huj limit we only want diagrams with at most two higher Landau level (n > 1) Green's 
functions. Green's functions meeting at a noise average must be in the same Landau level; while those meeting at a 
vertex can be in different Landau levels. So either the two n > 1 Green's functions meet each other at a noise average, 
or each meets three n = Green's functions at a vertex. Otherwise, there will be at least three higher Landau level 
Green's functions; see for example, Fig. 0(c). We have already considered those diagrams in which the two n > 1 
Green's functions meet each other at a noise average — they are the ones that yield the flux-flow result a xx (k) - 
they are also essentially two-point quantities and thus do not probe viscous effects such as shearing. 

Let us find an expression for those terms in which each n > 1 Green's functions meets a vertex. Recall the TDGL 
equation in symbolic form \I>i = G\$r)i — (5G\^2^i^2 (eq. ([L6|)). Of the four \&'s in the Kubo formula, let us replace 
the two that are in higher Landau levels with the second term on the right hand side of eq. (|16|). This substitution 
yields 



V xx (k Xl k z ) = ^*J m2 J d(r-r') J d(t-t')JJ 6 e^-*')^* 

(p lx + p 2 * x ) (p 3 , + p^S^G^Gs^nn^^n^c 



+G*G* 4 eO^tf*^**** 6 tf 3 >e + c.c. 



(57) 



l=2=(r,t) 
3=4=(r',t') 



In this symbolic notation, the numbers denote points in space-time, e.g. 5 — > (^,£5). This expression factors out 
the "bare" n > 1 Green's functions explicitly, leaving averages which involve LLL states only. The LLL average is 



a rather complicated eight-point object, but in principle, coupled with eqs. (49) and (p2|), it allows one to find the 
conductivity from simulations that use only LLL states. Because the Green's functions that have been factored out 
involve higher Landau level states, they are rather short-ranged in both time and space. When investigating a regime 
in which one expects the LLL object to become long-ranged, it is tempting to approximate the Green's function by 
some appropriately scaled delta function S(r — r')6(t — t 1 ) after the derivatives have acted. 

In the next section we will provide some explicit expressions for the diagrams in Figs. ^, and ^| for the two- 
dimensional case. But before leaving this section, let us make one point about the k z behavior. As previously 
mentioned the diagrams |i|(d) and ||(d) cancel at k z = 0, but note that in these diagrams there is no longer an n = 1 
channel through which we can pass the external momentum k z . Therefore, one might expect that one picks up the 
£ c length scale instead of £ c , however, this is not the case. Just as these diagrams cancel at k = so do the leading 
terms (i.e. the k z ^ terms) in the k z expansion. 



V. TWO-DIMENSIONAL RESULTS 



Many of the results above hold for films and even simplify in this case. The uniform conductivity is still given by 
the flux-flow formula in the LLL approximation. The diagrammatic structure of the perturbation series is exactly the 
same; and which diagrams are down and so forth also carries through; only the evaluation of the diagrams changes. 

To do these calculations we will of course need the 2D Green's function, which in the absence of an electric field is 

z J J n=0 n 

where E n = a + TiLjQ(n + 1/2) and L z is the film thickness. The Gaussian result for the conductivity is 
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o-xx(k) 



/2D 



£ 2 k 2 



h 2 LU%l 2 k 2 y 

8a 2 D 



Offc 4 ) 



(59) 



in the LLL limit. Other than a different numerical coefficient in front of the k 2 (S xyyx ), the result is similar to the 
3D result (eq. (f4(f)). Again we are applying the HF approximation so that a = a B + 2/3(|^o| 2 ) where 

k B T 



(|*0| 2 >2D 



which leads to the self-consistent equation 



OL2D = olh 



2tt a 2 D L z £ 2 ' 



P k B T 



7T £ 2 L Z a 2 D ' 

The expansion parameter in the two-dimensional perturbation series is 

/3 k B T 



X2D 



4tt PL z a 2 D ' 



allowing one to re-express eq. ( |6l|) as an — <52£>(1 — 4x2£>)@ 
Evaluating some of the diagrams in Figs. |^ and || gives 



a {FF) 

w XX 


= *£>(0) 


XX 




XX 


= 4?(o) 


^ XX 


= *£>(o) 


CT 6(a) 


= *£>(o) 


a m 

XX 


= *£>(o) 



7r 2 
< X 2D 



b 2D 



b 2D 



b 2D 



e -klP/2 ( 1 _ k 2J2 + l fe 2 fc 2 ^ 



- kli2/2 {i-k 2 y e + \k^ + k 2 x k 2 y ^ 



2 1)2 | 7,2 1)2 



-kit 



y 



h 2,2 pi , i.L2i,4 «6 
K y Kj_l -t- 



1 + /s 2 £ 2 + 2fc 2 £ 2 —kyk^i 4 + — k 2 k^£ 6 



r 2 
b 2D 



C -^J 2 ' 2 (kit 2 ~ k 2 y k 2 ± £ A + \k 2 yk\£^ 



e -3*j/72 ( _ fc 2 £ 2 + fc 2 £ 2 _ k 2 k 2ji + L k 2 k ^ 



(60) 



(61) 



(62) 



(63) 
(64) 
(65) 
(66) 
(67) 
(68) 



These calculations were done with the Green's function shown in bold in the figures in the n — 1 state. Recall that 
if we are interested in S xyyx we must redo the calculation with all n = states, and if we are interested in S xxxx we 
must also add to the expressions above those with n = 2 Green's functions. Here our interest is in showing explicitly 
the overall null contribution of these diagrams at k = — whether it be that the individual diagrams have zero 
contribution (e.g. diagrams 5(c), 6(a) and 6(b)) or that diagrams cancel (e.g. diagrams 4(d) and 5(d)) — and to 
suggest that the same does not apply at k ^ 0. 



VI. DISORDER 



Disorder introduces competing effects: it provides pinning centers which may couple to the viscous effects even at 
k = 0, but it also disrupts the growth in crystalline order which yields the viscous effects in the first place. With 
translational invariance destroyed, the flux-flow formula no longer applies. We will examine the effects of disorder at 
low order in perturbation theory within the Kubo formalism. 

We can in principle take the direct approach at k = using the Green's functions with the explicit electric field 
applied. However, it is more complicated than it was in the absence of disorder. We confirmed for the pure system 
that we could use solely n = states from the outset. So long as the electric field was explicitly applied, taking the 
ol/Tiloq — > limit in the beginning (i.e. using Green's functions with only n — states) or at the end yielded the same 
results. However, the same is not true in the presence of disorder. The time integrals here are more involved, and 
taking the a/hujQ — ^ limit in the beginning and at the end gives different results. With this complication, the direct 
approach is no better than using the Kubo formalism. 
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A. Point disorder 



In the experiments by Fendrich et a/0 the point defects induced by irradiation were seen to introduce an additional 
term to the ab plane conductivity which was not of the flux-flow form. They analyzed their resistivity data using 
the form p = (1/pf + l/pp)~~ 
latter had an activated form 

Let us consider the effect of uncorrelated disorder on the conductivity calculated from TDGL. We will model the 
point defects by adding a quenched random mass term to the free-energy functional 



1 where pf is the flux- flow resistivity and pp the resistivity due to point defects. The 
Pp = Po ex p{~U(T,H)/T}, and the field dependence of the "plastic energy" was 



d 3 i 



+ {a + a R (r))\V\ 2 + .. 



where the random variables a R (r) have zero average 



and ^-function correlations 



a R (r) = 



Wo 



ttfl(r) a R (r') = — S(t - r'). 
The associated disordered TDGL can be written as 

*(r,t) = f dv' [ dt'G(r,t;r',t') [r,(v', t') - a R (r')V(r', f) - /3|*(r', f')| 2 *(r', t')] 



or symbolically as 



*i = Gi. 2 r?2 - ajBGi, 2 f 2 - /3G X 0*2*2*2- 



(69) 



(70) 



(71) 



(72) 



(73) 



We have simply added the new term to eq. (|15|), hence G\^. remains the same, that is, the Green's function of the 
pure, linearized system. 

Iterating this equation with j3 = gives 



*i = Gi.irii - otrnGxpGi^z + amamGi^Gi^G^^riA + 0(a R ). 



(74) 



Diagrammatically, we will represent a R by a small square. There will be one arrow (Green's function) pointing into it 
and one out of it. The disorder averaging pairs up the a R s, we represent this feature by a dashed line connecting the 
two squares. Because of the spatial S function in the correlation of the disorder (|7l|), squares connected by a dashed 
line represent the same point in space but different points in time. Thus as far as spatial integrals are concerned, 
diagrams in this perturbation expansion are identical to those found expanding the nonlinear term. It is the time 
integrals that make it different. 

Figure [ij shows some diagrams that arise at the first order in an expansion in the disorder strength Wo- Note that 
if we were to draw the two points connected by the dashed line as one point, then we would end up with diagrams 
very much like those in Figure ||. (They would still differ by the number of circles.) In fact, the diagram in Fig. 0(d) 
gives zero contribution at k = for both a xx and a zz for the same reasons as that in Fig. ^(b). The diagram in Fig. 
M(a) is Green's function renormalizing, those in (b) and (c) noise renormalizing, and that in (d) vertex-renormalizing. 
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Certain basic facts about the calculations do not change in the presence of disorder. For instance, in the calculation 
oi<Jxx(k x , k y = 0, k z ), it remains true that there is a maximum of two n > 1 Green's functions in diagrams contributing 
to the a/fiujQ — > limit. Hence the diagram in Fig. 0(d) does not contribute to this order. Another invariant is that 
either the two n > 1 Green's functions meet each other at a noise average (as in Fig. 0(a) and (b)) or they each meet 
a vertex (as in Fig. 0(c)). The difference is now there are two types of vertices, the original quartic vertices and the 
new disorder induced vertices. From the case in which the n > 1 Green's functions meet at a noise average, one gets 
the disordered analog of the flux- flow formula with (|<lLd 2 ) replaced by its disordered counterpart. (This particular 
subset of diagrams yields the 1/pf in the Fendrich et al.t3 analysis.) From the case in which the two Green's functions 
meet quartic vertices, one gets the disordered analog of eq. (|57j). 

Then there are new cases, like the noise-renormalizing diagram shown in Fig. 0(c). Evaluating it yields 



i 



2 K z t c 



(75) 



Note that we have factored out the pure Gaussian result; it multiplies a numerical factor and zq, the dimensionless 
expansion parameteiH for point disorder 



Wn 



(76) 



We have also included the diagram's k z dependence. The c-axis length scale associated with this diagram, as well as 
those in Fig. (a) and (b), is i c . As in the pure case, we can send the external momentum k z through channel made 
entirely of n = 1 states. Recall that the points connected by the dashed line are the same point in space; consequently, 
in Fig. 0(c) the dashed line acts as a "short" allowing the external momentum k z to pass through a channel comprised 
solely of n > 1 states. 

By once again factoring out the higher Landau level Green's functions, this last calculation might be generalized to 



^ii (k x , k z ) — 



2 Wn 



16fc s Tm 2 



d(r - r') / d(t - t') / dr 



dt R 



dte. Q ik *( x ~ x ')+ ik z( z ~ z ') 



(Pi* + PL) (Pax + PL) |c*(r 2 , t; r 5 ,t 5 )G(r 3 ,t'; r 5 ,t 6 ) 
(*o(ri, t)%(r 5 , t 5 )*o(r 5 , t 6 )*5(r 4 , t')) c + c.c. I , 



(77) 



where we imagine the remaining LLL object as being fully renormalized. The disorder average was performed with 
the result that the order parameter fields v I'o( r 5, is) and v E , o( r 5i *5) are a ^ the same point in space but different points 
in time. Note that the LLL object here is a four-point quantity and thus probes viscous effects even at k = 0. The 
resummation suggested in eq. ([77]) is just one example of a new feature induced by the defects. 
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B. Correlated Disorder 



We model columnar defects lying parallel to the c axis by changing the correlation of the disorder from a 3D delta 
function to a 2D delta function 



a(r) a(r') 



Wl 



S(x - x') S(y - y'). 



(78) 



For sake of comparison we calculated the diagram in Fig. |?](c) for columnar defects, obtaining 



3zi 



(-2 + (1 + fc^/4) 1 / 2 + (l + ^ c 2 /4) 



3^2(1+^2/4)2 



(79) 



where we have factored the result as above and Z\ is the dimensionless expansion parameter for columnar defects 



z\ 



Wi 



(80) 



Note that z\ differs from zq by the absence of £ c in the denominator of zi, making z\ the larger of the two in the 
a^O limit. The c-axis length scale associated with the diagram in Fig. ^|(c) is £ c . Points connected by a dashed 
line are no longer at the same point in space. They have different values of z; therefore, the dashed line no longer 
serves as a "short" for external momentum k z and we can no longer send k z through an exclusively n > 1 channel. 
The switch from the £ c length scale for point defects to £ c for columnar defects in the evaluation of theses diagrams 
is in accord with the suggestion that correlated defects align pancake vortices enhancing their integrity. 
We model planar defects by changing the disorder correlations to 



a(r) a(r') 



W 2x 



8{x — x'); 



a(r) a {v') = ^L 8{y-y'), 



(81) 



where the former models planar defects parallel to the yz plane and the latter those parallel to the xz plane. When 
the electric field is in the x direction, the vortices move in the y direction. One might expect different results for 
vortices moving parallel to or perpendicular to the defects. Evaluation of the diagram in Fig. |?](c) is then 



<7^)(fc z , )=a^)(0) 



r 7( c ; 



(0) 



8-y/2 

i\/^Z2y 

8V2 



(-2 + (l + A:^ c 2 /4) 1 / 2 + (l + fc^ c 2 /4) 

3^(1 + ^2/4)2 
(-2 +(1 + £^/4)l/2 + (l + fc2 4 2/ 4 ) 



3^2(1 + ^2/4)2 



where 



Z2i = 



W 2 r 



1671-te 2 



(82) 



(83) 



(with i — x or y) are the dimensionless expansion parameters for planar defects. For this particular diagram we 
find no difference between the xz and yz planar defects. The planar defect expansion parameter z 2 i has the same 
a dependence and so the same temperature dependence as z%. On the other hand, zi has an additional £ in the 
denominator so the magnetic field dependence of these expansion parameters is different. 



VII. SUMMARY 



We have investigated the conductivity in the vortex-liquid regime via perturbation theory. Using an approach 
in which an electric field is applied explicitly we showed that within the lowest Landau level (LLL) approximation 
the uniform ab plane conductivity, a xx (0), is proportional to (l^ol 2 ) with a temperature- independent proportionality 
constant, to all orders in perturbation theory. This result elevates the derivation of the flux-flow formula to non- 
zero temperatures. We identified the subset of diagrams that yields the same result calculated within the Kubo 
formalism and conjectured that the remaining diagrams must cancel. We verified this cancellation to second order in 
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perturbation theory and also demonstrated that the cancellation does not extend to the nonlocal conductivity a xx (k) , 
allowing viscous effects to enter. 

When it comes to the nonlocal conductivity and the issue of characteristic lengths, the present work has certainly 
raised more questions than it has answered. For instance, the lowest-order (Gaussian) diagram of the transverse ab- 
plane conductivity already revealed a long, temperature-dependent length; while all of the diagrams for the longitudinal 
a6-planc conductivity individually had short, temperature-independent length scales. If the longitudinal length scale 
is to grow long, as is suggested by its identification with shearing effects of a growing crystal, it will be through the 
combined effect of many diagrams. The nature of the diagrams that may lead to such an effect have been identified, 
but what remains unclear is whether the longitudinal and transverse length scales are independent or related. 

There is a similar question regarding the c-axis. Every diagram contributing to the LLL approximation of a zz 
has £ c , the temperature-dependent correlation length, as its length scale. On the other hand, at each order in the 
perturbative expansion of a XXl the c-axis length scale is £ c , a temperature-independent magnetic length. If a xx {z — z') 
is to become substantially nonlocal as T is lowered, it will be through the combined effect of various orders. If the two 
length scales are long, it is unclear whether they will turn out to be essentially the same or distinct. In addition, there 
is the question of whether or not the a6-plane and c-axis length scales grow independently. Finally, these matters 
must be readdressed with regard to the nonlocal resistivity which may have its own distinct length scales. 

We also examined the effect of disorder on the afr-plane conductivity. While there was a subset of diagrams that lead 
to the disordered analog of the flux-flow formula, there were other contributions that coupled to the viscous effects 
even at k = 0. We also showed that diagrams contributing to the LLL approximation had as their c-axis length scale 
£ c in the presence of point defects and £ c in the presence of columnar and planar defects. 

Missing from this work is any consideration of what vector-potential fluctuations might do. Including them would 
require a new kind or vertex, one involving a $, a "J* and an A^, where A^ is the fluctuating part of the vector 
potential. The arguments leading to the flux-flow result may no longer hold when these new vertices are added. Also 
lacking are the effects of Maxwell's equations which are needed to form a complete set of equations. At low order 
these lead to back-flow effects and at higher orders may produce lattice effects. There are many challenging aspects 
of this problem to be resolved. 
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